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$_i^ ' The origin of the atmosphere of the largest moon of Saturn, Titan, is 

J-»N ' poorly understood and its chemistry is rather complicated. Ground-based 

j^ ^ millimeter/sub-millimeter heterodyne spectroscopy resolves line shapes suffi- 

ciently to determine information in Titan's atmospheric composition (on ver- 
tical profiles and isotopic ratios) . We test the capabilities of the Swedish Het- 
erodyne Facility Instrument (SHFI), Receiver APEX-1, together with the At- 
^ acama Pathfinder Experiment APEX 12-m telescope for Titan's atmospheric 

0> , observations. In particular we present sub-millimeter observations of the CO(2- 

Cn . 1) and HCN(3-2) lines of the Titan stratosphere with APEX, and with SHFI 

•/j ' taken during the Science Verification (SV) instrument phase on March and 

C^ ' June 2008. With the help of appropriate radiative transfer calculations we in- 

vestigate the possibility to constrain the chemical concentrations and optimize 
^~^ • the performance of the APEX-1 instrument for inferring vertical profiles of 

^^v ' molecular components of the atmosphere of Titan. 

This study attempts to contribute to constrain radiative transfer and retrieval 
algorithms for planetary atmospheres, and to give hints to the current and 
future ground and space-based data acquisition leading to a more thorough 
1 J understanding of the chemical composition of Titan. 

r\ , 
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sponding to the observing proposal number 081.F-9812(A). 

'Currently at the National Institute of Information and Communications Technology, 
Japan 



May 21, 2010 2:32 WSPC - Proceedings Trim Size: 9in x 6in titan-rsh-reresubmission-vS 



1. Introduction 

Titan's atmosphere exhibits a complex photochemistry. The origin of car- 
bon monoxide (CO) is not well understood (whether photochemical or pri- 
mordial). Hydrogen Cyanide (HCN), the most abundant nitrile in Titan, 
is a key intermediate in production of more complex hydrocarbons and 
organic molecules. HCN in the Titan atmosphere has been discovered by 
the infrared observations of Voyager 1^, and CO has been detected by the 
ground-based near-infrared observations^. Following these detections, the 
molecular concentrations of CO and HCN in Titan's atmosphere have been 
determined from infrared, millimeter, and sub-millimeter observations as 
well as from modeling (see Tables 2 and 3 in section 4 for details). CO seems 
to be uniformly mixed in Titan's atmosphere up to high altitudes. HCN 
abundances, however, display a steeper profile with ambiguous enrichment 
values (Tables). Already several observational data confirm that Titan's 
atmospheric composition is indeed seasonally^ '^, and spatially dependent^. 
Full behavior of the seasonal characteristics of the spatial distribution is 
not constrained yet due to the limitation in the temporal coverage of the 
previous observations, and to the use of different instruments. New sub- 
millimeter observations arc required not only to provide new abundance 
constraints and shed more light on the rate of these seasonal variations, 
for example, but for support, complement and cross-calibration the ESA's 
Herschel Space Observatory mission. One of the goals of the Key Program 
of Herschel, entitled Water and Related Chemistry in the Solar System^, 
is to understand water inventory in the Titan's atmosphere as well as dis- 
tributions of other hydrocarbons and nitriles. During the guaranteed time, 
line surveys on Titan at the frequency range of 500 and 5000 GHz are go- 
ing to be carried out using two low-to-medium resolution spectrometers 
(Photodetector Array Camera and Spectrometer (PACS) and Spectral and 
Photometric Imaging Receiver (SPIRE)). Our model calculations of the 
synthetic spectra of Titan expected to be observed with Herschel show 
that several CO and HCN lines will be detected with PACS and SPIRE'^. 
Because of their low-to-mcdium spectral resolutions, SPIRE and PACS are 
not capable of measuring the shapes of CO and HCN lines, which prevents 
us from determining precisely the vertical profile of CO and HCN mixing 
ratios. Therefore, high spectral resolution observations of CO and HCN 
with ground-based sub-millimeter telescopes would significantly improve 
the accuracy of retrieving the CO and HCN profiles. 

The Atacama Pathfinder Experiment (APEX) 12-m telescope is oper- 
ational since 2005®, and has been already used to the planetary science 
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(monitoring mesospheric winds of Venus^), among numerous galactic and 
extra-galactic objects. In spring 2008, the new APEX Swedish Heterodyne 
Facility Instrument (SHFI) has been commissioned on the APEX 12-m tele- 
scopc"'^'^'^^, which consists of four single-pixel receivers (APEX-1, APEX-2, 
APEX-3, and APEX-T2) mounted in a single cryostat located in the Nas- 
mythA cabin of APEX. In this report we present a summary of our ob- 
servations of CO and HCN performed with APEX-1 during the Science 
Verification (SV) phase with the aims to demonstrate the capabilities of 
APEX and SHFI for Titan's atmospheric observations, and to investigate 
the possible retrieval of CO and HCN abundances. We also examine the 
possible detection of the millimeter/sub-millimeter rotational transition of 
C2H2 (already observed in the mid-infrared spectra by the Composite In- 
frared Spectrometer, CIRS, on Cassini^^), and the search for one new com- 
pound, HC5N(83-82) which the possible presence has been suggested by 
experimental results simulating Titan's atmosphere i3-i5^ 

The observations arc described in detail in section 2 while the radiative 
transfer modeling and the data analysis arc discussed in sections 3 and 
4, respectively. A brief discussion is given in section 5 and a summary in 
section 6. 

2. Observations 

Titan's spectra were acquired with the APEX 12-m telescope located on 
Chajnantor (Chile) during 21 March, and 17-18, 23, and 27 June 2008 dur- 
ing SV. We used the APEX-1 receiver, operating at 211-270 GHz, which 
employs a superconductor-insulator-superconductor (SIS) mixer and be- 
haves as a single sideband receiver (SSB). For the backend, we used the 
Fast Fourier Transform Spectrometer (EFTS) with a channel separation 
of 122kHz and a bandwidth of 1 GHz. Tablet summarizes the measured 
transitions as well as the observing days and the precipitable water vapour 
(P WV) . Titan was observed near the western or eastern elongations at sep- 
aration angles from Saturn greater than 120 ". Pointing and focusing of the 
telescope were regularly checked scanning across Saturn in azimuth and 
in elevation (APEX has a pointing accuracy of 2" r.m.s. over sky). The 
apparent diameter of Titan was around 0.8". 

Because a bug in the low-level Python ephemerides software not pre- 
viously detected. Titan was not tracked properly at the beginning of our 
observations on 21 March 2008. The bug was successfully removed by the 
APEX staff. These observations in particular have helped to improve the 
control software. Accounting for possible frequency errors in the line center 
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Table 1. Species observed on Titan 



Species 



Frequency (GHz) Date 2008 PWV (mm) 



HCN(3-2) 
HC5N(83-82) 
C0(2-l) 
HCN(3-2) 
C2H2 J=47 vb= 



J=46 1/4=1 



265.886 


21 March 




220.937 


17-18 June 


3.5-4.0 


230.538 


23 June 


0.5 


265.886 


23 June 


0.5 


254.521 


27 June 


0.2 



(of a few MHz around the respective line centers), and considering the short 
integration time, we could not detect any evident signature of HC5N(83-82) 
and C2H2 in the recorded spectra averages of integration times of 36 and 
91 min. (r.m.s. of 0.035 and 0.03 K), respectively. Lower r.m.s. noise levels 
of around 0.01 K (which could be achieved with integration times of around 
8 and 13 h, respectively) would have facilitated the possible detections of 
these species. Therefore we concentrate here on the C0(2-l) and HCN(3-2) 
observations. Because of the high oversubscription of observing proposals 
and the uncertainties about the tracking at the beginning of the campaign, 
the initially proposed integration time of around 2 hours for C0(2-l) and 
2.65 hours for HCN(3-2), which corresponds to a r.m.s noise level of 0.017 K 
and 0.025 K respectively, was not available. Instead, C0(2-l) and HCN(3-2) 
were integrated for around 13 and 19 min, respectively. The initial reduc- 
tion of each spectrum was performed using the CLASS software package of 
the Grenoble Astrophysics Group**. 

Fig. 1 presents the recorded Titan's spectra of the C0(2-l) and HCN(3- 
2) lines. The resolution has been downgraded to about 4 MHz in order to 
increase the signal-to-noise ratio. 

The conversion from the observed antenna temperature T* to the total 
flirx density from the Titan's atmosphere Sv,tot [Jy] is performed by S^^tot = 
24 Ax T* ■ rjf/rja, where rjf and rja are the forward and aperture efficiencies, 
respectively. In this study, we adopt the values considering the efficiencies 
listed on the APEX web site*' and of the instrument paper^: rja are 39.9 and 
40.2 for 230 and 265 GHz, respectively, and ?/y is 0.95 for both frequencies. 
These values are expected to have a 10% uncertainty. 

By fitting the synthetic spectra to the observed one, we intend to infer 
the abundance of CO and HCN in Titan's stratosphere. A precise compar- 
ison of the observations with synthetic spectra must cope with two limita- 



''URL:http://www.iram.fr/IRAMFR/GILDAS 
'^URL:http://www. apex-telescopc.org/tolescope/efficiency/ 
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Fig. 1. Titan's spectra observed with the APEX 12-m Telescope on 23 June 2008. 
Left: CO(2-l) line rebinned over 32 resolution elements (4 MHz). Right: HCN(3-2) line 
rebinned over 32 resolution elements (4 MHz). 



tions: (1) the linc-to-continuum ratios of these CO and HCN spectra were 
not able to be confirmed from the observed data. 1 GHz bandwidth is too 
narrow for constraining the continuum IcveL We have examined the possi- 
bility of using the emission at the frequency of ±500 MHz offset from the 
line center as the pseudo-continuum. Wc found, however, that the intensi- 
ties of the CO (2-1) line at the ±500 MHz significantly depend on the CO 
abundance at a certain range of altitude, which also affects the intensity at 
the line center. Therefore, we conclude that using such fluxes as the pseudo- 
continuum is not a suitable method for determining the CO abundance in 
this study. This is also the case for the HCN (3-2) line. (2) Baseline ripples 
and uncorrected larger scale baseline characteristics may affect the shape 
of the observed spectra. One can obtain a spectral line shape symmetric 
with respect to the line center by folding the spectrum. However, such a 
spectra folding may introduce an artificial error in the line shape. Therefore 
we prefer to keep the spectra with their original baseline distortions. 



3. Radiative transfer modeling 

The synthetic spectra are generated by a multi-layered line-by-line radiative 
transfer model considering the spherical geometry of Titan's atmosphere. 
It consists of 120 layers that span the 0-600 km interval with a resolution 
of 5 km. In general, the atmospheric model for Titan is explained fully in 
Rengel et al. (2009)^. Briefly, Titan's thermal and pressure profiles based on 
the in-situ measurement by the Huygcns probc^^ are adopted for the opac- 
ity calculation of the CO and HCN gases at each altitude level. The mixing 
ratio profiles of CO and HCN are set as the parameters to be retrieved 
in this study. The detail of modeling vertical profiles of each molecule 
is described later. For CO, the assumed N2-broadening line-width here is 
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that of Semmoud-Monnanteuil & Colmont (1987)^'', and for HCN, the value 
corresponding to the same transition at the infrared vibrational-rotational 
band^^ as there is no exact measurement for this pure rotational line. The 
collision induced absorption coefficients of N2; CH4, and H2 mixtures, which 
dominate the millimeter and sub-millimeter continuum opacity of Titan's 
atmosphere, are included by using the formulation of Courtin (1988)^^ and 
Borysow and Tang (1993)2°. 

4. Data analysis 

4.1. The CO abundance 

Recent observations suggest that CO in Titan is uniformly mixed through- 
out Titan's atmosphere^^"^^. Such a profile can be supported if its long 
chemical lifetime in the oxygen-poor Titan atmosphere is considered ('--'10"^ 
years^^). In this study, we retrieve CO abundance by assuming a vertically 
constant profile. It should be noted that some of the previous observations 
have shown the possibility of existence of the discontinuity in the CO mix- 
ing ratios between the troposphere and the stratosphere (e.g., ~30ppm 
in the troposphere^^ and ~60ppm in the stratosphere^^). However, such 
discontinuity cannot be constrained by our observations, as the CO (2-1) 
measurement reported here does not have any sensitivity to the altitude 
below 60km (Fig. 2 left panel). 

A direct comparison of synthetic and observed spectra is made and the 
best fitting model with the least x^ is determined. In Fig. 2 right panel, 
the best-fit CO abundance shows a solution of BOlg^ ppm. Although this 
solution agrees with the previous stratospheric CO measurements of Hi- 
dayat et al. (1998)^^, Lellouch et. al (2003)^5, and Baines et al. (2006)^2 
(Table 2), we have to allow a larger error on our result by considering the 
uncertainty in the fiux scaling and the baseline problems. To examine the 
effect of the uncleaned baseline of the observed spectrum, we performed 
the fitting analysis by constraining the data only on each side of the wings. 
Using the red-shift wing {v < vq) yielded the best fitting CO mixing ratio 
of 26^8^ ppm, while using the blue-shift wing {v > vq) resulted 36^^^ ppm. 
Taking into account the 10% uncertainty in the fiux calibration, the ac- 
ceptable range of CO mixing ratio becomes as 16-51 ppm. This wide range 
covers most of the past works. One way to avoid such a large ambiguity in 
the data analysis (i.e. to avoid the effect of the uncertainty in the flux cali- 
bration) is to calibrate the observed spectrum by scaling its continuum level 
to that of the forward model calculation. Our forward model calculations 
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suggest that the intensity at the frequency offset larger than 
be used as the continuum level. 
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Fig. 2. Left: Weighting functions with respect to CO mixing ratio for selected frequen- 
cies near the CO(2-l) line center. Note that the finite spectral resolution (4 MHz) of 
the measurement is taken into account for the calculation of these weighting functions. 
Right: Observed and synthetic CO(2-l) spectra at 230.54 GHz (grey spectrum and line 
ones, respectively). The synthetic spectra are calculated assuming vertically uniform CO 
mixing ratios of 10, 30, 50, and 70ppm (from bottom to top) respectively. The small 
upper panel shows a zoom of the spectrum for a frequency range of ±60 MHz . 



4.2. The HCN vertical profile 

We retrieve the mixing ratio profile of HCN by performing forward- 
and inversion calculations iteratively. For solving the inversion problem, 
we employ the optimal estimation method (OEM) (e.g., see details in 
Rodgers (1976)^^). The Levenberg-Marquardt scheme is used in the iter- 
ation. We set the logarithm of HCN mixing ratio as the retrieving param- 
eter. The a priori profile of HCN mixing ratio, which is required in the 
OEM approach to regularize the ill-posed inversion problem, is based on 
the previous ground-based work by Marten et al. (2002)'^'^. We assumed 
the a priori error of 2.3 in the logarithm scale. During the HCN retrieval, 
the temperature profile is fixed to our nominal temperature profile (the one 
based on the Huygens measurements). 

Fig. 3 shows the best-fit synthetic spectrum (solid line) and the corre- 
sponding HCN profile. The averaging kernels of the retrieved HCN mixing 
ratios indicate that this retrieval has a significant sensitivity at the alti- 
tude around 90 km though we can still retrieve information at the higher 
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Table 2. Observations of CO in Titan atmosphere 



Altitude 


Mixing ratio [ppm] 


Wavelength* [/J-m] 


Facility/Instrument 


Reference 


Troposphere 


40+100 


1.57 


KPNO-4m/Fouricr Transform interferometer 


2 


Stratosphere 


30-180 


35314.89 


VLT-70m/27 antennas 


28 


Stratosphere 


60±40 


2602.59 


OVRO/2 elements interferometers 


29 


Stratosphere 


2tl 


2602.59 


IRAM-30m/3-mm SIS receiver 


30 


Stratosphere 


50±10 


2602.59 


OVRO/6 10.4 m diam. antennas 


31 


Troposphere 


lOi^o 


4.8 


UK IR Telcscope/CGS4 spectrograph 


32 


Stratosphere 


27± 5 


2602.59, 1304.34, 900.9 


IRAM-30m/SIS heterodyne receivers 


27 


Stratosphere 


52±6 


1304.34 


OVRO/6 antennas 


33 


Troposphere 


32±10 


4.75-4.85 


VLT 8-m/ISAAC 


25 


Stratosphere 


51±4 


869.56 


SMA/5 and 6 antennas 


21 


153-350 km 


60 


4.5-4.85 


VLT 8-m/ISAAC 


26 


Tropo-Stratosphere 


45±15 


4.64 


Cassini/CIRS 


34 


Stratosphere 


32±15 


4.64 


Cassini/VIMS 


22 


Stratosphere 


47±8 


333.3-166.66 


Cassini/CIRS 


35 


Stratosphere 


301^5 


1301.29 


APEX/APEX-1 receiver 


This article 



*In spite that each wavelength region is defined by standard spccifical units , we use here 
fira to facilitate the comparisons. 



altitudes up to 240 km. Our retrieved HCN mixing ratios are larger than 
those of Marten et al (2002) by a factor of --30-80% at the altitudes of 100- 
200 km. The synthetic spectrum calculated with the a priori HCN profile, 
i.e., the result of Marten et al. (2002), is also shown for comparison (Fig. 3). 
A large discrepancy between the observed and the synthetic spectra, and 
the a priori profile is seen mostly at the wings of the spectrum. This results 
in the retrieved larger HCN abundances at the altitudes of 90-150 km. Note 
that our result showing lower HCN abundance values at around 300 km is 
out of the altitude range where our observations are sensitive. 

In the above mentioned retrieval, we only took the system tempera- 
ture of the APEX-1 instrument into account for the statistical error of the 
measurement. In fact, as we discussed previous in section 4.1, the observed 
spectrum contains the 10% uncertainty in the flux calibration. If a 10% 
larger scaling factor is applied in the flux calibration, we cannot derive any 
realistic HCN profiles which can reproduce the observed spectrum. In case 
of using a 10% smaller scaling factor, the HCN abundances for the best-fit 
spectrum are larger than those of Marten et al. (2002) at the altitude of 
100-150 km (which is attributed by the wings of the spectrum) and then 
becomes as small as ~30% of Marten et al. (2002) at higher altitudes of 
150-200 km (Fig. 4). 

The two retrievals (Figs. 3 and 4) demonstrate the significant impact of 
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the flux calibration uncertainty on the retrieved HCN profile. Since we can- 
not constrain the line-to-continuum ratio of the observed spectrum, we are 
not able to conclude the HCN abundance at altitude of 150-200 km which 
deviates within a wide range from 30 to 150% of the a priori guess. For 
the lower altitudes such as 100-150 km, the wing of our observed spectrum 
requires a higher opacity of HCN than that of Marten et al. (2002) even if 
we consider the uncertainty in the flux calibration. 
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Fig. 3. Left: Observed and synthetic HCN(3-2) spectra at 265.89 GHz: Solid lino shows 
the best-fit spectrum, and dash line the spectrum calculated with the a priori HCN 
profile. Grey line is the observed data. The upper panel figure represents the zoom to 
the frequency range of -|-/-100MHz. Right: The retrieved HCN vertical profile (solid 
line) under the nominal temperature profile without considering the error in the fiux 
calibration. The error bar indicates the 1-tr limit of the retrieval. The dotted profile 
represents the a priori values where the gray shaded region corresponds to their errors. 
The ratio of the retrieved profile with respect to the a priori profile is also shown. On 
the extrem right panel, the averaging kernels of the retrieved HCN mixing ratios are 
shown. 
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Fig. 4. See caption Fig. 3, but here the measurement data is scaled by a factor of 0.9. 
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Table 3. HCN mixing ratios in Titan's atmosphere 



Altitude 


Mixing ratio [ppm] 


Wavelength* [/im] 


Facility/Instrument 


Reference 


Stratosphere 


3.0 X 10-'' 


3386.0 


IRAM-30m/SIS receiver 


38 


Stratosphere 


1.6 X 10-'^ 


14.025 


Voyagcr/IRIS 


39 


Stratosphere 


(0.75-52) X 10-'' 


3386.0 


IRAM 30-m/SIS receiver 


40 


Stratosphere 


4.7x 10-**-1.5 X lO"'' 


14.025 


Voyager/ IRIS 


41 


400 km 


~ 2 X 10-5 




Model prediction 


42 


Stratosphere 


(0.5-4) X lO-'^ 


14.025 


IRAM 30-m/SIS receiver 


43 


400 km 


1 X lO"'' 




Model prediction 


44 


700 km 


1 X 10-5 




Model prediction 


44 


400 km 


1 X lO-f' 


3386.0 


IRAM-30m/SIS receiver 


45 


400 km 


1 X 10^5 




Model prediction 


46 


700 km 


1 X 10-4 




Model prediction 


46 


Stratosphere 


3.0 X lO-'' 


14.025 


ISO/SWS 


47 


~600 km 


7 X 10-3 


3 


Keck II/NIRSPEC 


46 


83 km 


3x 10-5 


1692.43 


SMA/5 and 6 antennas 


21 


300 km 


3 X 10-5 


1692.43 


SMA/5 and 6 antennas 


21 


500 km 


0.4 X 10-«* 


14.04 


Cassini/CIRS 


48 



*In spite that each wavelength region is defined by standard specifical units , we use here 
fim to facilitate the comparisons. 



5. Discussion 

How does one explain the behavior of the retrieved HCN vertical profile ob- 
tained here? the relative poor signal-to-noise of the observations obtained 
here results in a poor fit to the HCN line wings. Furthermore, comparing 
the results obtained with previous results derived from disk-averaged infor- 
mation, interferometer measurements, and space-based observations sum- 
marized in Tables must be considered cautiously. In other words, direct 
comparison of the results from different instruments requires us to assume 
that all instrument-related offsets have been accounted for. 

Certainly, Titan's atmospheric temperature and composition are not 
in a steady state. Does the HNC abundance retrieved here perhaps indi- 
cate that the composition of Titan's atmosphere has changed during the 
last years? Fig. 7 in Teanby et al. 2009^ shows that the HCN abundance 
measured by Cassini (at southern latitudes) has not significantly changed 
between 2006 and 2008. Then it seems to be that possible detections of 
time variations (if any) between different observations could be maybe due 
to observing more of the northern hemisphere than of the southern one as 
Titan's season changes. As we derived disk-averaged information we need 
to keep in mind, however, that our vertical distributions are mostly repre- 
sentative of the equatorial region of Titan since the measured flux density 
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Fig. 5. Comparison of recent observational and modeling works of HCN abundances. 
The small numbers correspond to the indices of the reference on Table 3. The shaded 
region represents the results here obtained, considering the l-cr limit of the retrieval 



spectra are more heavily weighted toward equatorial latitudes. In case that 
vertical distributions would be representative of the northern latitudes, the 
profile analysis would become rather complicated. There is evidence of a 
vortex, which acts to separate enriched from unenriched air, at 25 — 55 ° N 
encircling the North Pole 5.34,49^ 

Further higher signal-to-noise HCN observations at different times are 
necessary in order to open the possibility to detect seasonal and even spa- 
tial variations in HCN, and confirm or not an enrichment layer in Titan's 
stratosphere. 



. Summary 

We report the first observations obtained with the SFHI APEX-1 instru- 
ment on APEX on a planetary/satellite atmosphere taken during SV. It 
consists of the spectra of C0(2-I) and HCN(3-2) in Titan's atmosphere. 
The observations reported here improved the control software of the 
APEX telescope, so that now it is possible to track planets with the 
APEX-1 receiver. 

We investigate the CO and HCN composition of Titan's stratosphere. 
Our CO mixing ratio estimation is consistent with other authors. We 
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retrieved the HCN vertical profile, which is inconsistent with previous 
analyses. This requires further investigation by future HCN observations 
with higher signal-to-noise ratios. The search of nitriles and CO appears 
favorable in the sub-millimeter range explored with the APEX telescope 
and with the APEX-1 receiver. 

This work represents a first step in exploring the capabilities of APEX 
and SHFI for Titan's atmospheric observations. New simultaneous obser- 
vations of HCN and its isotopes at higher frequencies with APEX-2 and 
APEX-T2, for example, would improve the sensitivity to the abundance 
retrievals, and potentially could support, complement, and cross-calibrate 
the ESA's Herschel Space Observatory mission. 
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